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ABSTRACT 

We report on high-resolution spatial and timing results for binary millisecond pulsar PSR J0218+4232 
obtained with the Chandra HRCT and HRC-S in imaging mode. The sub-arcsecond resolution image 
of the HRC-I (0.08-10 keV) showed that the X-ray emission from PSR J0218+4232 is consistent with 
that of a point source excluding the presence of a compact nebula with a size of about 14" for which 
we had indications in ROSAT HRI data. The presence of a DC component is confirmed. This X-ray 
DC component has a softer spectrum than the pulsed emission and can be explained by emission from 
a heated polar cap. With our HRC-S observation we obtained a 0.08-10 keV pulse profile with high 
statistics showing the well-known double peaked morphology in more detail than before: The two pulses 
have broad wings and pulsed emission appears only to be absent in a narrow phase window of width 

0.1. The absolute timing accuracy of ~ 200/is makes it possible to compare for the first time in 
absolute phase the X-ray pulse profile with the highly structured radio profile and the high-energy 7-ray 
profile (0.1-1 GeV). The two X-ray pulses are aligned within absolute timing uncertainties with two of 
the three radio pulses. Furthermore, the two 7-ray pulses are aligned with the two non-thermal X-ray 
pulses, corresponding to a probability for a random occurrence of 4.9cr, strengthening the credibility of 
the earlier reported first detection of pulsed high-energy 7-ray emission from a (this) millisecond pulsar. 

Subject headings: pulsars: individual (PSR J0218+4232 ), X-rays: stars 



1. INTRODUCTION 

PSR J0218+4232 is a 2.3 ms pulsar in a two day orbit 
around a low mass (~ 0.2 Mq) white dwarf companion 
(Navarro et al. 1995; van Kerkwijk 1997). Pulsed X-ray 
emission with a Crab-like double pulse profile has been 
reported from ROSAT 0.1-2.4 keV data (Kuiper et al. 
1998a) and BeppoSAX MECS 1.6-10 keV data (Mineo et 
al. 2000). The pulsed spectrum as measured by the MECS 
appeared to be remarkably hard with a power-law photon 
index 0.61 ± 0.32, harder than measured for any other ra- 
dio pulsar. Furthermore, Kuiper et al. (2000) report the 
detection with EGRET of pulsed high-energy (0.1-1 GeV) 
7-ray emission from this ms pulsar. Their argument is 
based on three lines of evidence: (1) the 0.1-1 GeV data 
show a 3.5<7 pulsation at the radio period; (2) the 7-ray 
light curve resembles the one seen in hard X-rays, namely 
a phase separation of ~ 0.45 between two pulses/maxima; 
(3) the spatial analysis shows that the position of the 
EGRET source 3EG J0222+4253 moves from the position 
of the nearby BL Lac 3C 66A towards the pulsar position 
with decreasing gamma-ray energy (for energies between 
100 and 300 MeV all source counts could be attributed to 
PSR J0218+4232). They also showed that the two 7-ray 
pulses/maxima appeared to be aligned in absolute phase 



with two of the three radio pulses detected at 610 MHz. 
Unfortunately, the timing accuracies of the ROSAT and 
BeppoSAX observations were insufficient to construct X- 
ray profiles in absolute phase. None of the current mod- 
els for pulsed X-ray and 7-ray emission from radio pul- 
sars offers a consistent explanation for the above summa- 
rized high-energy results on PSR J0218+4232 (Kuiper et 
al. 2000). 

PSR J0218+4232 is also remarkable in that it is the only 
Crab-like ms pulsar with a large DC (unpulsed) fraction 
of 63 ± 13% in the ROSAT band below 2.4 keV (Kuiper 
et al. 1998a), as well as a large DC fraction of ~ 50% 
in radio, systematically over the range 100-1400 MHz 
(Navarro et al. 1995). The DC components as measured 
in the ROSAT and radio observations could be explained 
by emission from a compact nebula with diameter ~ 14", 
but in both cases the indications were at the limit of the 
imaging capabilities. Assuming that the radio DC compo- 
nent is compact, combined with the measured very broad 
and structured radio pulse profile, Navarro et al. (1995) 
suggested that the magnetic field of PSR J0218+4232 is 
almost aligned with the rotation axis, the observer viewing 
the system under a small angle with respect to the rotation 
axis. Stairs et al. (1999) measured the magnetic inclination 
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Fig. 1.— (left) Chandra 0.08-10 keV HRC-I image of a 10" X 10" region centered on the radio pulsar position of PSR J0218+4232 . The 
radio position is marked with a '+' sign. The angular distance between the radio pulsar position and the X-ray centroid is ~ 0"6, consistent 
with the Chandra localization accuracy, (right) Radial distribution of HRC-I (Level 1) events using the best fit maximum likelihood X-ray 
position as centre. Superposed as dotted line is the radial profile from the PSF determined from the 2D fit procedure, which is consistent 
with the Chandra response to a point source. The dashed line indicates the background level derived from counts in the range 10 - 25 arcsec 
from the centre. 



angle analyzing radio polarization profiles. Their rotation 
vector model fits indicate that the magnetic inclination an- 
gle is indeed consistent with 0° (8° ±11°). Unfortunately, 
in their fits the line-of-sight inclination is unconstrained. 
If the DC component in X-rays is also compact, for the 
suggested geometry of a nearly aligned rotator and a small 
viewing angle, it can originate in the pulsar magnetosphere 
as well as from a heated polar cap of the neutron star. 

The objectives of our Chandra observations were: 1) To 
establish the spatial extent of the X-ray DC component, 
compact or extended; 2) To construct an X-ray pulse pro- 
file which can be compared in absolute phase with radio 
profiles and the 3.5cr EGRET high-energy (0.1-1 GeV) 7- 
ray profile. 

2. OBSERVATIONS 

PSR J0218+4232 was observed with the HRC-I and 
HRC-S instruments on the Chandra X-ray observatory 
(CXO) during two observations of <~ 75 ks duration each. 
The HRCs are multi-channel plate (MCP) detectors sen- 
sitive to X-rays in the 0.08-10 keV energy range with no 
spectral information. 

The first observation with the HRC-I camera at the fo- 
cal plane of the telescope mirror took place on 1999 De- 
cember 22 for an effective exposure of 74.11 ks. Unfor- 
tunately, this observation suffered from a non-recoverable 
timing ( "wiring" ) problem assigning incorrectly the event 
trigger time to that of the next trigger (Tcnnant et al. 
2001). On-board event screening prevents recovering the 
true event times of all the triggers in the telemetry stream 
by simple back-shifting. The timing accuracy in this case is 
determined by the total event rate, typically 300-400 s _1 , 
and is much worse than the 16/zs intrinsic timing resolu- 
tion. This degradation prevents the construction of high- 
resolution pulse profiles in case of millisecond pulsars. 



In a re-observation of PSR J0218+4232 on 2000 Octo- 
ber 5 for an effective exposure of 73.21 ks with the HRC-S 
in imaging mode (using only the central MCP segment) 
all event triggers could be put in the telemetry stream, so 
that back-shifting recovers the intrinsic timing accuracy. 

3. SPATIAL ANALYSIS WITH HRC-I 

PSR J0218+4232 has clearly been detected near the cen- 
tre of the 30'x30' field of view in the 74.11 ks HRC-I obser- 
vation. Within 14'5 from the pulsar 51 X-ray sources have 
been detected, some of them have counterparts in our pre- 
vious ROSAT HRI (0.1-2.4 keV) and BeppoSAX MECS 
(1.6-10 keV) observations (Kuiper et al. 1998a; Mineo et 
al. 2000). The central ~ 7' x 7' part of the X-ray image con- 
taining PSR J0218+4232 has also been observed at opti- 
cal wavelengths by the Keck telescope (van Kerkwijk 2001) 
and it turns out that all central sources except one have op- 
tical counterparts. Relative astrometry demonstrates that 
the Chandra positions are accurate to within 0'.'5(1), bet- 
ter than the celestial localization accuracy determination 
requirement of 1". Using the standard Level 1 event file 
provided by the CXC pipeline analysis (this file contains 
all HRC triggers with the position information corrected 
for instrumental (de-gap) and aspect (dither) effects), a 
zoom-in at the pulsar location (Fig. 1 left) shows that the 
X-ray centroid has a ~ 0'.'6 offset from the radio-pulsar 
position of 2 h 18 m 6?351, 42 32'17'.'45 (epoch J2000) consis- 
tent with our previous findings. It is also evident from the 
image shown in Fig. 1 (left) that the source events are 
concentrated within a circle with a radius of <~ 2" from 
the X-ray centroid. A more quantitative estimate involves 
maximum likelihood fitting of the measured 2D spatial 
distribution with a model composed of a Gaussian and a 
generalized Lorentzian with free shape parameters, nor- 
malization and centroid position on top of an (assumed) 
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Fig. 2.— Pulse-profile of PSR J0218+4232 as measured 
by the HRC-S in imaging mode. Two cycles are shown for 
clarity. Typical error bars are shown at phases 0.35 and 0.7. 
The DC-level (±lcr) is indicated by the hatched area, and 
the background level from the imaging analysis by the dashed 
line. 



flat background, also with free scale factor. The optimum 
X-ray centroid position is 2 h 18 m 6?305(l),42°32T7'.'23(2) 
(epoch J2000) while the radial profile of the best model 
fit is compatible with the PSF of the HRMA/HRC-I com- 
bination (95% of the source counts are within 2" from 
the X-ray centroid). Fig. 1 (right) shows the best model 
profile (dotted line) superposed on the measured radial 
profile using the optimum X-ray centroid position as cen- 
tre. Thus, we have no evidence for extended emission near 
PSR J0218+4232 at - 1" scales (diameter), rejecting the 
indication for a compact nebula found in our analysis of 
ROSAT HRI data (Kuiper et al. 1998a). 

The total number of counts assigned to the pulsar from 
the maximum likelihood fit is 870 ± 30, which translates 
in a count rate of (1.175 ± 0.041) x 10~ 2 . Assuming ab- 
sorption in a column with density Nh = 5 • 10 20 cm ~ 2 
(Verbunt et al. 1996) and a photon power-law index of 
0.94 (Mineo et al. 2000) the count rate converts to a 2-10 
keV flux of (5.2 ± 0.2 ± 2 ;J) x 10~ 13 erg cm^s- 1 . The 
first error represents the statistical error and the second 
the systematic error due to uncertainties in the absolute 
sensitivity of the HRC-I and in the spectral model. The 
value is consistent with the BeppoSAX MECS value of 
(4.38 ± 0.48 ± 0.44) x 10~ 13 erg cm^s" 1 (the first error 
represents the statistical uncertainty and the second error 
the systematic uncertainty of about 10%). 

4. TIMING ANALYSIS WITH HRC-S 

In a spatial analysis using the HRC-S Level 1 event 
file the X-ray centroid of the X-ray counterpart of PSR 
J0218+4232 was found to have an offset of 1'.'9 from the ra- 
dio position of the pulsar, larger than the celestial location 
accuracy determination requirement of 1"; a systematic 
difference noticed earlier by e.g. Wang et al. (2001). The 
number of counts assigned to the source, applying a simi- 
lar likelihood procedure as in the case of the HRC-I image 
data, is 755 ± 28. Using the optimized background scale 
factor and source profile, the optimum radius for event ex- 
traction (S/N optimal) from the best fit centroid turns out 



to be l'/5. 

The first step in the timing analysis is to correct the 
assigned event times by back-shifting, recovering the in- 
trinsic relative time resolution of 16/xs. Next, the event 
extraction radius was set to the optimum radius of 1'.'5. 
The final step is the determination of the arrival times at 
the solar system barycentre using the orbital information 
of Chandra and the position of PSR J0218+4232. Folding 
the barycentered arrival times with the spin and binary pa- 
rameters from an updated ephemeris for PSR J0218+4232 
revealed the well-known doubled peaked profile at high 
statistics (Fig. 2). The deviation from a flat distribution 
is 15.2cr according to a Z| - test and the peak separa- 
tion is 0.475 ± 0.015 consistent with previous estimates 
(Kuiper et al. 1998a; Mineo et al. 2000). The peak widths 
A<j> (FWHM), fitting the profile in terms of two asym- 
metric Lorentzians atop a flat background, are: Peak-1 
(at 4> ~ 0.2) 0.116 ± 0.019 (267 ± 44/xs) and Peak-2 (at 
<j> ~ 0.7) 0.109 ± 0.025 (251 ± 58yus). Fitting two Gaus- 
sians plus background yielded similar results. 

In Fig. 2 is indicated an estimate for the unpulscd 
(DC) level (±lc) using a bootstrap method outlined by 
Swanepoel et al. (1996). Furthermore, the background 
level estimated in the maximum likelihood imaging anal- 
ysis is also shown, giving a number of background counts 
within the l'/5 extraction radius of (46.9 ± 0.4). The 
pointlike (see Sect. 3) DC (unpulscd) component, sig- 
nificantly detected with ROSAT (0.1-2.4 keV; 4.8cr) by 
Kuiper et al. (1998a) and weakly with BeppoSAX (1.6- 
4.0 keV; <~ 2a) by Mineo et al. (2000) is clearly seen in 
this 0.08-10 keV profile. The pulsed fraction Tp defined 
as Afp I {Afp + Afvc) , where Afp specifies the number of 
pulsed source counts and AT-dc the number of DC source 
counts, turns out to be Tp = 0.64 ± 0.06 as measured by 
CXO HRC-S over the entire 0.08-10 keV energy range. 

Finally, our CXO HRC-S 0.08-10 keV X-ray pulse profile 
can directly be compared in absolute phase with our previ- 
ously derived radio and 7-ray profiles (Kuiper et al. 2000). 
Firstly, because we generated an ephemeris valid for all 
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Fig. 3. — Multi-wavelength pulse profiles of PSR 
J0218+4232 in absolute phase, a) Radio pulse profile at 610 
MHz. b) Chandra HRC-S X-ray pulse profile (0.08-10 keV; 
timing accuracy 200^is or 0.085 in phase) and c) EGRET 7- 
ray pulse profile (0.1-1 GeV; timing accuracy 100/^s or 0.043 
in phase). The absolute timing accuracies of the X- and 7-ray 
profiles are shown as horizontal bars centered on phase 0.62. 
Indicated as dotted lines are the positions of the 2 pulses 
in the 610 MHz radio profile which coincide with the high- 
energy pulses. Typical ±lcr error bars arc indicated in the X- 
and 7-ray profiles. Note that the background level in the 7- 
ray profile as determined from an imaging analysis is at 22.2 
(Kuiper et al. 2000). 
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used observation periods and applied a common reference 
epoch (see Table 1), and secondly, the absolute timing ac- 
curacy of CXO is within 200^s (Tennant et al. 2001) which 
corresponds to a phase inaccuracy smaller than 0.085 given 
the pulse period of 2.3 ms. Fig. 3 shows then for the first 
time the comparison of a PSR J0218+4232 X-ray profile in 
absolute phase with the radio and 7-ray profiles. Within 
the uncertainties of each of the measurements the X-ray 
pulses are aligned with two of the three radio pulses at 610 
MHz and with the 7-ray pulses. 

5. CHANDRA HRC-S X-RAY PULSE PROFILE VS. EGRET 
HIGH-ENERGY 7-RAY PROFILE 

The alignment in absolute phase of the very-hard X-ray 
pulses with the low-significance 7-ray pulses is particularly 
interesting. The calculated 3.5ct significance of the 7-ray 
profile was derived for unbinned data giving the statistical 
evidence {Z\ - test) for deviations from a flat distribution 
anywhere in pulsar phase. When we calculate the signif- 
icance of detection of the two 7-ray pulses at the known 
absolute phases of the X-ray pulses (e.g. by determining 
the pulsed excess counts in the narrow pulse phase win- 
dows 0.10 - 0.30 and 0.55 - 0.70, considering the remaining 
phase windows (65%) as background; we found 54.8 ± 13.7 
pulsed excess counts) then we derive for this single trial 
an increased total significance of 4. Oct. 

In fact, we performed simulations similar to what we 
did earlier for PSR B1951+32 (Kuiper et al. 1998b): We 
simulated 10 6 flat phase distributions with the same num- 
ber of counts (n = 308) as contained in the EGRET 0.1-1 



GeV phase histogram shown in Fig. 3c. Then we ap- 
plied the Z\ - test and verified first that the distribution 
of the Z\ simulations behaved as expected for a random 
distribution (with 10 6 trials we could calibrate well up 
to ~ 4.5ct). Next, we used again the two narrow pulse 
phase windows 0.10 - 0.20 and 0.55 - 0.70 to determine the 
pulsed excess counts. Fig. 4 shows in a two-dimensional 
distribution all simulations with a Z\ - significance above 
3ct versus the measured pulsed excess counts in the pre- 
defined pulse windows. As can be seen, only for one sim- 
ulation out of 10 6 a significance above 3.5ct and pulsed 
excess counts above 55 was reached. This corresponds to 
a statistical significance of ~ 4.9ct. Therefore, we take the 
alignment in absolute phase of the non-thermal X-ray and 
7-ray pulses as supporting evidence for our first detection 
of high-energy 7-rays from a millisecond pulsar. 

6. SPECTRAL ANALYSIS 

In the spatial and timing analyses we derived the main 
results from our Chandra data, which contain very limited 
spectral information. Essentially, we measured the inte- 
gral count rate over the entire 0.08-10 keV energy window. 
Therefore, we cannot add significant new information to 
the well-defined pulsed spectrum as measured by the Bcp- 
poSAX MECS instrument (Mineo ct al. 2000). However, 
no meaningful information has been published sofar on 
the spectral shape of the DC component (crucial to decide 
on the origin and production mechanism) , while the com- 
bined ROSAT HRI, BeppoSAX MECS and Chandra data 
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Fig. 4. — Scatter plot of pulsed excess counts vs. Z\ for 
10 6 background simulations. The hatched band indicates the 
parameter space of interest, namely, to detect by chance a sig- 
nal with a modulation signifcance above 3.5(7 and more than 
55 pulsed excess counts in the pulse phase windows 0.1 — 0.3 
and 0.55 — 0.7 (representative for the 0.1-1 GeV 7-ray profile). 
We found just one occurrence in 10 6 trials, corresponding to 
a statistical significance of ~ 4.9o\ 
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can give a first estimate: The ROSAT HRI provides us a 
DC count rate for 0.1-2.4 kcV, the MECS for 1.6-4 keV 
and 4-10 kcV, and Chandra for 0.08-10 keV. 

Assuming a spectral model and absorbing column Nh 
and using the instrument sensitivity curves over the dif- 
ferent energy windows, the HRC-S DC count rate can be 
used to predict the ROSAT HRI and BcppoSAX MECS 
DC count rates. 

For a black body model the measured ROSAT 
HRI/HRC-S rates demand a temperature kT <, 0.9 kcV; 
the measured BeppoSAX MECS/HRC-S rates 0.5 <, kT 
<5 1.05 (ler limits). Therefore, a thermal model with a 
temperature 0.5 ^5 kT ^ 0.9 (90% confidence interval: 
0.33-1.1 keV) is consistent with all the measured DC count 
rates. Assuming a power-law model, the DC rates of the 
three instruments are consistent with power-law indices 
in the range 1.3-1.85 (la limits; 90% confidence interval 
0.95-2.45), significantly softer than for the pulsed emis- 
sion (0.61 ± 0.32) as measured by the MECS (Mineo et al. 
2000). The total statistics of the data are insufficient to 
discriminate between the black body and the power-law 
spectral shapes. 

7. DISCUSSION & CONCLUSIONS 

Our first objective of the Chandra observations was to 
determine the spatial extent of the X-ray DC component 
and to find a likely explanation for its origin. As is shown 
in Fig. 1, no evidence for extended emission near PSR 
J0218+4232 at ~ 1" scales (diameter) is seen, rejecting 
the indication for a compact nebula found in our anal- 
ysis of ROSAT HRI data (Kuiper et al. 1998a). The 
combined ROSAT, BeppoSAX MECS and Chandra data 
led us to conclude that the DC spectrum is significantly 
softer than the pulsed spectrum, the latter being clearly 
of non-thermal origin. If we adopt the geometry proposed 
by Navarro et al. (1995), a nearly aligned rotator and a 
small viewing angle, the soft DC X-ray component can 
be explained as thermal emission from the polar cap of 
the neutron star which stays visible to the observer all 
the time. For old millisecond pulsars the thermal emission 
can only originate from reheating of the polar cap area by 



back-flowing accelerated particles. In fact, for a ms-pulsar 
the polar-cap half angle 9 pc (— arcsin(-\/(i? n;s • 2n/p ■ c))) 
extends over a relatively large angle, e.g. 17?5 for PSR 
J0218+4232 . This means that the spin-axis is well within 
the (magnetic) polar-cap region, for the measured an- 
gle between the spin- and magnetic axes amounts ^10° 
(8° ± 11°) (Stairs et al. 1999). Therefore, a small observer 
viewing angle to the system (angle spin-axis and line of 
sight) can keep the observer looking at the heated polar 
region continuously. 

The Navarro et al. (1995) radio observations had a VLA 
beam size of 16", also allowing an interpretation of the ra- 
dio DC component (reported DC fraction ~ 50%) as a 
compact nebula up to the size of the VLA beam. No new 
observations have been reported with a smaller beam size, 
but Stairs et al. (1999) and Kuzmin & Losovsky (2001) 
revisited the pulsed emission. We fitted the total radio 
spectrum of Navarro et al. (1995) (five data points) with 
a power-law shape, index 2.57 ± 0.07, and the new total 
pulsed radio spectrum (three data points from Navarro 
et al. (1995), two from Stairs et al. (1999) and one from 
Kuzmin & Losovsky (2001)), index 2.58 ±0.15. The spec- 
tra appeared identical in shape and the updated DC frac- 
tion is (13 ± 9)%, significantly lower than the earlier esti- 
mate. This leads us to the conclusion that the radio DC 
component is most likely also compact and has the same 
magnctosphcric origin as the pulsed emission, which has a 
very structured and remarkably broad profile and is indeed 
practically never "off" . 

Our second objective was to obtain with Chandra a 
more significant X-ray profile than measured sofar, to 
study the profile structure in more detail in the X-ray band 
below 10 keV, and to have for the first time absolute tim- 
ing, allowing multiwavelengths phase comparisons. The 
profile in Fig. 2 with a significance of 15.2a is indeed 
much more significant than the MECS profile (6.8c). Now 
we can really see that the two pulses have broad wings and 
that the profile reaches the DC level only in a very narrow 
phase interval around phase 0.35 (phase extent & 0.1). 
In fact, Fig. 3 shows that the radio and the X-ray profiles 
both reach a minimum level for approximately the same 
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absolute phase region. Furthermore, Fig. 3 shows that the 
X-ray pulses are aligned in absolute phase with two of the 
three radio pulses and the two 7-ray pulses within the tim- 
ing uncertainties of the different measurements. 

PSR J0218+4232 is the first and only millisecond pulsar 
for which we reported evidence for detection of high-energy 
gamma-ray emission up to 1 GeV (Kuiper et al. 2000). We 
take the alignment in absolute phase of the non-thermal 
X-ray and 7-ray pulses as important supporting evidence 
for our first detection of high-energy 7-rays from a mil- 
lisecond pulsar. In Sect. 5 we showed that the probability 
that a random timing signal reaches a 3.5u modulation 
significance and has its 55 pulsed excess counts in phase 
with the two X-ray pulses amounts <~ 4.9o\ 

Wallace et al. (2000) reported the detection of short- 
term variability in the high-energy 7-ray emission from 
3EG J0222+4253 (3C66A/PSR J0218+4232) in a sys- 
tematic search for all 170 unidentified sources in the 3rd 
EGRET Catalog (Hartman et al. 1999). They noted that 
if this flaring is due to PSR J0218+4232 : "It would be an 
unusual source in two ways: it would be the only millisec- 
ond pulsar seen by EGRET and the only pulsar to show 
strong flaring" . We argue, however, that in their search 
the reported evidence for flaring is not significant: For 
each of the 170 unidentified sources Wallace et al. (2000) 
produced light curves per viewing period (VP) with 2 day 
flux values each time they were in the field-of-view during 
a VP. For 3EG J0222+4253 (3C66A/PSR J0218+4232) 
they found the most significant evidence for variability, 
namely in VP 15 with a variability index V=2.6, corre- 
sponding to a ~ 3cr significance for a random detection in 
a single trial assuming Gaussian statistics (although Pois- 
son statistics apply) . This variability was due to one single 
high 2 day flux value. However, this source was viewed in 
4 observations (no indication for variability in the other 
VPs) which makes the probability (4 x 0.002512) 1% or 
2.58(7 to find V=2.6 in 1 out of 4 observations. Further- 
more, Wallace et al. (2000) analysed 144 VPs with dura- 
tion larger than 3 days producing for the 170 unidentified 
sources few- hundred light curves, making the probability 
to find once V=2.6 (or one single high 2 day flux value) few 
times unity! They correctly noted that their Monte Carlo 
probabilities give misleading values (too optimistic) as the 
used averages having no uncertainties assigned. Also in 
the latter case the very large number of "trials" has been 
ignored. Since there was no a priori reason to select 3EG 
J0222+4253 for a single trial, we donot regard the indica- 
tion for short-term variability of this source significant. 

The X-ray results on the DC emission and the radio 
results on the DC and pulsed emission suggest for PSR 
J0218+4232 an emission scenario of a nearly aligned ro- 
tator and a small viewing angle. The latter angle can 
still be as large as ~ 20° to explain the X-ray DC emis- 
sion if this originates from a polar-cap heated by parti- 



cle bombardment. For explaining the pulse profile with 
two hard-spectrum X-ray/7-ray peaks it is more critical 
to know this viewing angle to get a handle on the ge- 
ometry. For both competing classes of models, polar cap 
models (PC; e.g. Daugherty & Harding (1994, 1996)) and 
outer gap models (OG; see e.g. Cheng et al. (1986a, b); 
Ho (1989)) production of hard X-ray/7-ray emission in 
the magnetospheres of millisecond pulsars has been pre- 
dicted (e.g. Bhattacharya & Srinivasan (1991); Sturner 
& Dermer (1994)). Recently, millisecond pulsars among 
them PSR J0218+4232 were considered for PC models by 
Dyks & Rudak (1999); Bulik & Rudak (1999); Bulik et al. 
(2000); Zhang & Harding (2000). For PSR J0218+4232 
they did not succeed in reproducing the measured high- 
energy spectral shape. More recently, Dyks & Rudak 
(2002) succeeded in reproducing the high-energy spectrum 
of PSR J0218+4232, but required non-orthodox assump- 
tions about the electron energy distribution or emission 
altitude as well as off-beam viewing geometry. Wozna 
et al. (2002) reproduced the double-peak profile with the 
measured phase separation for a small inclination angle of 
8° and a viewing angle of 29° (angles approximately con- 
sistent with our findings). However, for these angles the 
spectrum reached its maximum luminosity for too high 
energies around 100 GeV. The fact that in the case of 
PC models a double-peak hard X-ray/7-ray profile can be 
modelled for a nearly aligned rotator with a ~ 20° viewing 
angle follows already from Daugherty & Harding (1996), 
who showed that the 7-ray beam produced at the polar cap 
rim is relatively wide, amounting ~ 26° (1.5 x 9 pc ) for PSR 
J0218+4232. For emission produced at higher altitudes in 
the magnetosphere (as is proposed in the more recent ver- 
sions of PC models), the beam will become even wider. 
In all these discussions the real viewing angle in combi- 
nation with the inclination angle determines whether the 
measured double peak profile can be obtained. Also for 
OG models a double peak profile is in principle possible, 
depending on the actual viewing angle and which part of 
the outer gap (in altitude) is visible to the observer. For 
OG models the hard X-ray/7-ray cone might become too 
broad for the case of PSR J0218+4232 as a nearly aligned 
rotator if the hard X-ray/7-ray production takes place only 
close to the light cylinder. However, Hirotani et al. (2002) 
argued recently that the 7-ray production is not limited to 
regions in the outer gap above the null-charge surface, but 
can also originate from altitudes closer to the neutron star. 
More detailed model calculations are required to solve the 
present uncertainties in the overall interpretation of the 
data. 

This underlines the importance of a new attempt to de- 
termine the geometry from radio polarization data. Then 
detailed model calculations starting from this geometry 
can attempt to reproduce all the remarkable timing and 
spectral results obtained for this millisecond pulsar. 
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Table 1 

Ephemeris of PSR J0218+4232 based on radio timing measurements spread over a 6.5 year period 



Parameter 


Value 


Right Ascension (J2000) 


02 h 18 m 6*351 


Declination (J2000) 


42° 32' 17"45 


Epoch validity start/end (MJD) 


49092 - 51462 


Frequency 


430.4610670731491 Hz 


Frequency derivative 


-1.43397 x 10~ 14 Hz s" 1 


Epoch of the period (MJD) 


50277.000000025 


Orbital period 


175292.302218 s 


a • sin i 


1.9844301 (lt-s) 


Eccentricity 





Longitude of periastron 





Time of ascending node (MJD) 


50276.61839862 


RMS of timing solution 


14.4 milliperiods; 33.4 /is 



